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■ ABSTRACT 
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Context. Seismic modelling of the f} Cephei stars promises major advances in our understanding of the physics of early B-type stars on (or 
close to) the main sequence. However, a precise knowledge of their physical parameters and metallicity is a prerequisite for correct mode 
] identification and inferences regarding their internal structure. 
"■---«. Aims. Here we present the results of a detailed NLTE abundance study of nine prime targets for theoretical modelling: y Peg, 5 Cet, v Eri, /3 
^ ' CMa, f 1 CMa, V836 Cen, V2052 Oph, p Cep and DD (12) Lac (hereafter 12 Lac). The following chemical elements are considered: He, C, N, 
O, Mg, Al, Si, S and Fe. 

Q ■ Methods. Our curve-of-growth abundance analysis is based on a large number of time-resolved, high-resolution optical spectra covering in 

' most cases the entire oscillation cycle of the stars. 

CZ) . Results. Nitrogen is found to be enhanced by up to 0.6 dex in four stars, three of which have severe constraints on their equatorial rotational 

velocity, Q.R, from seismic or line-profile variation studies: /? Cep (Q.R~26 km s '), V2052 Oph (Qft~56 km s '), S Cet (Q.R < 28 km s ') and 
* * i I i 

^> £ CMa (ClRsini S 10 km s ). The existence of core-processed material at the surface of such largely unevolved, slowly-rotating objects is 
• »"H . not predicted by current evolutionary models including rotation. We draw attention to the fact that three stars in this subsample have a detected 
magnetic field and briefly discuss recent theoretical work pointing to the occurrence of diffusion effects in j3 Cephei stars possibly capable of 
altering the nitrogen surface abundance. On the other hand, the abundances of all the other chemical elements considered are, within the errors, 
indistinguishable from the values found for OB dwarfs in the solar neighbourhood. Despite the mild nitrogen excess observed in some objects, 
we thus find no evidence for a significantly higher photospheric metal content in the studied ji Cephei stars compared to non-pulsating B-type 
stars of similar characteristics. 
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1. Introduction 

The class of the /3 Cephei pulsating variables is defined as 
B0-B3 V-III stars excited by both (low-order) pressure and 
gravity modes. The concommitant existence of these two types 
of modes, along with their relatively simple internal struc- 
ture, offers the prospect of probing their deep interiors via de- 
tailed seismic modelling. This will allow stringent constraints 
on some fundamental parameters of early B-type dwarfs and 

Send offprint requests to: T. Morel, e-mail: 
thierry@ster . kuleuven . be. 

* Table IA. II is only available in electronic form at the CDS via 
anonymousftptocdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat?J/A+A/7??/??? 
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subgiants, such as the extent of convective core overshooting 
or the rotation law in the radiative envelope, to be set in the 
near future. Although forthcoming or space missions already 
in operation (e.g. MOST, COROT) are expected to lead to dra- 
matic advances in this field, very intensive ground-based obser- 
vations have already demonstrated the potential of such tech- 
niques (Dupret et al. l2004> . 

The unstable modes in [} Cephei stars are driven by the k- 
mechanism and arise from an opacity bump at T ~ 2 x 10 5 
K in their interior. It immediately follows that the incidence 
of pulsations is predicted to be largely controlled by the metal 
content, with the instability domains in the Hertzprung-Russell 
(HR) diagram shrinking or even completely vanishing in the 
low metallicity regimes (Pamyatnykh 119991 Deng & Xiong 
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2001 1. However, recent pieces of evidence suggest that our un- 
derstanding of the excitation mechanisms in these stars is likely 
to be still incomplete. First, the number of ft Cephei candi- 
dates in the LMC is much higher than theoretically predicted 
(Kolaczkowski et al. l2004> . Second, current pulsation models 
have increasing difficulties in accounting for all the frequen- 
cies detected in some key objects as more and more intensive 
observational campaigns are undertaken (e.g. 12 Lac: Handler 
et al. 120061 v Eri: Pamyatnykh et al. 120041 Ausseloos et al. 
2004 1. Non-standard stellar models including gravitational set- 
tling and radiative levitation might need to be invoked in such 
cases, as first suggested by Pamyatnykh et al. (2004i. 

These inconsistencies between theory and observations will 
be better appraised when the physical parameters (especially 
r e ff) and chemical composition (i.e. stellar opacities) of these 
objects are precisely known. Noteworthy attempts to deter- 
mine the effective temperature and global metallicity of several 
ft Cephei stars from International Ultraviolet Explorer (IUE) 
data have recently been presented (Niemczura & Daszyriska- 
Daszkiewicz 2005 1, but more robust estimates are likely to be 
obtained from a detailed analysis of high-quality optical data. 
Furthermore, a knowledge of the abundances of the individual 
chemical species is preferable in the context of theoretical mod- 
elling, as early-type stars may present significant departures 
from a scaled solar mixture. In particular, rotationally-induced 
mixing can lead to the dredge up of some core-processed CNO 
material in fast-rotating, main-sequence OB stars (e.g. Meynet 
& Maeder 120001 Proffitt & Quigley 2001 1. Such deviations 
from the standard solar abundance pattern might be routinely 
incorporated in the near future in any evolution and pulsation 
codes, and their impact on the pulsational properties explored 
(see Ausseloos 2005 1. 

Notwithstanding the ease with which high-quality spectro- 
scopic data can be gathered for such bright objects, surpris- 
ingly little attention has been paid in recent years to derive the 
stellar parameters and abundances of ft Cephei stars from high- 
resolution optical spectra, with the only dedicated studies go- 
ing back to the early 70s (Watson fTWTl fT972l Peters fl973l . 
It is the goal of this project to remedy this situation by tak- 
ing advantage of the dramatic improvement in data quality, at- 
mospheric models and non-LTE (NLTE) modelling techniques. 
Here we present a pilot study of nine archetypical ft Cephei 
stars, describing in detail the methodology used to derive the 
atmospheric parameters and elemental abundances in a self- 
consistent way. The next papers in this series will be devoted 
to a sample of Slowly Pulsating B stars (SPBs) and to the OB- 
type stars within reach of the asteroseismology programme of 
the COROT mission (Baglin l20"03l . 

2. Observational material 

Most of the selected stars have been the subject of very inten- 
sive (multisite) campaigns over the last few years, both in pho- 
tometric and spectroscopic modes. As such, they are amongst 
the ft Cephei stars (and indeed early-type stars) with the high- 
est number of pulsation modes identified and are thus ideally 
suited for further theoretical modelling (e.g. v Eri: Ausseloos 
et al. 120041 V836 Cen: Aerts et al. 120031 . Table □ summarizes 



Table 1. Basic pulsation properties of our targets. A^(v): num- 
ber of independent pulsation frequencies detected. Data from 
stated papers and references therein. 



Star 


Dominant pulsation 


N(v) 


Most recent 




mode 




reference 


y Peg 


Radial 


4 


Chapellier et al. 1 2006 1 


<5Cet 


Radial 


4 


Aerts et al. 1 2006 1 


v Eri 


Radial 


14 


Jerzykiewicz et al. 1 2005 l 


/3CMa 


Non-radial 


3 


Desmet et al. <2006t 


CMa 


Radial 


1 


Saesen et al. <2006> 


V836 Cen 


Non-radial 


6 


Aerts et al. J2004b» 


V2052 Oph 


Radial 


2 


Neiner et al. <5lJU3l 


/3Cep 


Radial 


5 


Telting et al. fT997l 


12 Lac 


Non-radial 


11 


Handler et al. J20051 



the basic pulsation properties of our targets. Most of the stars 
with a dominant radial mode do not possess large-amplitude 
non-radial modes, except v Eri. On the other hand, the radial 
mode is excited in the three stars with a dominant non-radial 
mode, but with a lower amplitude. 

A description of the spectroscopic data used is given in 
Table [2] Further details on the data acquisition and reduc- 
tion procedures can be found in Aerts et al. d2004al v Eri), 
Mazumdar et al. d2006l ft CMa), Saesen et al. d2006l £ CMa), 
Aerts et al. J2004bl V836 Cen), Neiner et al., in prep. (V2052 
Oph) and Desmet et al., in prep. (12 Lac). The data for 6 
Cet were obtained at ESO during the 2002 July 17-25 pe- 
riod and were reduced with a dedicated reduction software 
for the CORALIE spectrograph (see Baranne et al. 1996 1. 
The spectra of y Peg and ft Cep have been retrieved from 
the ELODIE archives (see Moultaka et al. I2004> . Aerts et al. 
( 2004b I presented a preliminary spectral synthesis analysis of 
V836 Cen based on a spectrum reduced with the FEROS reduc- 
tion pipeline. The raw data have been completely re-reduced 
here with standard IRAF 1 tasks, resulting in a much higher data 
quality. 

Our analysis for each star is based on a mean spectrum gen- 
erally created by co-adding a large number of individual expo- 
sures (up to 579; Table |3. These spectra were put in the lab- 
oratory rest frame prior to forming the weighted (by the S/N 
ratio) average of the whole time series, as large radial velocity 
variations are often observed. The orders of the mean spectrum 
were then merged using IRAF tasks. Sections of the merged 
spectrum were subsequently continuum normalized by fitting 
the line-free regions with low-order cubic spline polynomials. 
Each individual echelle order was separately normalized for the 
data acquired with the CS2 spectrograph at McDonald obser- 
vatory, as the merging procedure proved unsatisfactory in that 
case. 

The very high signal-to-noise ratio often attained for the 
combined spectrum enables us to confidently measure weak 
diagnostic lines. Figure ^ illustrates this point for some Fein 
features and allows one to get a sense of the overall data qual- 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 2. Basic description of the spectroscopic data. Spectral types from the SIMBAD database (revised values are given in 
Table[7}. AA: spectral range covered in A, R: mean resolving power of the spectrograph, yV: number of exposures combined, S /N: 
typical signal-to-noise ratio at about 4500 A in the individual exposures (estimated from photon statistics). 





y Peg 


5Cet 


v Eri 


/3CMa 


f 1 CMa" 
Max EWs Min EWs 


V836 Cen 


V2052 Oph 


/JCep 


12 Lac 


HD number 


HD 886 


HD 16582 


HD 29248 


HD 44743 


HD 46328 


HD 129 929 


HD 163 472 


HD 205 021 


HD 214993 


Spectral type 


B2IV 


B2IV 


B2III 


B 1 II/III 


Bl III 


B3 V 


B2 IV-V 


B2 Illevar 


B2III 


Telescope 


1.9-m OHP 


1 .2-m Euler 


1.2-m Euler 


1.2-m Euler 


1 .2-m Euler 


2. 2-m Euler 


2.7-m McDonald 


1.9-mOHP 


2.7-m McDonald 


Instrument 


ELODIE 


CORAL IE 


CORALIE 


CORALIE 


CORALIE 


FEROS 


CS2 


ELODIE 


CS2 


AA 


3895-6820 


3875-6820 


3875-6820 


3875-6820 


3875-6820 


3555-9215 


3630-10 275 


3895-6820 


3630-10275 


R 


50 000 


50 000 


50000 


50000 


50 000 


48 000 


60000 


50 000 


60000 


N 


47 


4 


579 


449 


38 41 


1 


105 


28 


31 


S/N 


215 


150 


250 


250 


180 200 


150 


230 


215 


285 



a Min and Max EWs correspond to minimum and maximum EWs of the Si m lines, respectively. 




- Fe III 528a. 3 

5280 5290 5300 5310 

Wavelength (A) 

Fig. 1. Mean, normalized spectrum of /3 CMa in the spectral 
range 5279-53 14 A. Some weak Fe m lines are indicated. 

ity. More importantly, combining such a large number of time- 
resolved spectra also allows us to minimize the impact of the 
strong line-profile variations arising from pulsations (see Aerts 
et al. 1 19941 for some examples) and the related changes in the 
equivalent widths (EWs), which are, for instance, typically of 
the order of ?>-\% for the Sim lines in f3 Cephei stars (De 
Ridder et al. 2002 1. All the fundamental quantities we derive 
(e.g. T e fi) can thus be regarded in most cases as representative 
of the mean values averaged over the pulsation cycles. This is 
one of the most important aspects of the present work with re- 
spect to previous studies based on snapshot spectra (e.g. Gies 
& Lambert [T9"9"2l hereafter GL). In the case of f 1 CMa, the ex- 
cellent time sampling and relatively simple nature of the pulsa- 
tions (only a large-amplitude radial mode is excited) prompted 
us to examine the phase-related changes in the atmospheric pa- 
rameters. To this end, we separately analyzed two mean spectra 
created for two phase intervals of width A0=O. 1 corresponding 
to minimum and maximum EWs of the Si ra lines (see Saesen 

et ai. renna . 

3. Method of analysis 

3.1. Model atmosphere calculations and spectral 
synthesis 

We made use of the latest versions of the NLTE line forma- 
tion codes DETAIL and SURFACE originally developed by 
Butler ( 1984 1 and Giddings d!981> . along with plane -parallel, 
fully line-blanketed LTE Kurucz atmospheric models (Kurucz 
I1993> . Grids with solar metallicity were used. Slight deviations 



from the solar mixture for some metals (e.g. N; see below) do 
not have an appreciable effect on the atmospheric structure. We 
adopted models with He/H=0.089 by number in all cases, in ac- 
cordance with the derived helium abundances. The only excep- 
tion was V2052 Oph, as it displays some evidence for a helium 
enrichment (see Sect.|4]and Neiner et al. l2003> . 

We used Kurucz Opacity Distribution Functions (ODFs) 
with a fine frequency mesh (1212 points). To limit compu- 
tational time, we adopted a constant microturbulent velocity, 
£mod=8 km s _1 , in the statistical equilibrium calculations with 
DETAIL. This value is in most cases consistent with the micro- 
turbulence derived from the line analysis. Test calculations for 
/3 CMa show that the exact choice of this parameter has little 
impact on the final results, with the mean abundance differ- 
ences being less than 0.05 dex for all elements when adopting 
fmod=2or8kms _1 . 

The following chemical elements are considered: He, C, N, 
O, Mg, Al, Si, S and Fe. The number of levels for the various 
ionic species is provided in Table|3] along with the original lit- 
erature sources where basic information on the model atoms 
can be found. Note that the atomic data used (e.g. photoion- 
ization cross sections) have been regularly updated following 
advances in theoretical calculations. The Fe in model atom will 
be described in detail in a forthcoming paper (Butler et al., in 
prep.). Here we note that the bulk of the atomic data, energy 
levels, oscillator strengths and damping parameters are taken 
from Kurucz ( 1993 1. While the Fe n and Fe iv atoms are rudi- 
mentary, the Fern model includes all terms calculated to lie 
below the ionization threshold. 

The line atomic data (e.g. oscillator strengths, radiative 
damping constants) are taken from the NIST 2 or VALD (Kupka 
et al. 1999 1 databases. The NLTE populations are computed as- 
suming LS coupling. Our line list includes most of the diagnos- 
tic lines generally used in abundance analyses of early B-type 
stars, but great care has been taken to only retain unblended 
features for the relevant temperature ranges. 3 Our final line list 
is made up of about 180 spectral lines. 



2 Available online at: http : //physlab2 .nist . gov/PhysRefData/ 
ASD/index.html. 

3 We made extensive use of the spectral atlases for main-sequence 
B stars available at: http://www.lsw.uni-heidelberg.de/ 
cgi-bin/websynspec . cgi (see Gummersbach et al. 1998 1. 
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Table 3. Number of levels for the various ions having their pop- 
ulations explicitely treated in NLTE with DETAIL. The ground 
states of C rv, N rv, O iv, Mg m, Al iv, Si v, S iv and Fe v are also 



considered. 


Ion 


Number 
of levels 


Reference 


Hi 


10 


Husfeld etal. (1989i 


He i 


27 


Husfeld et al. Il989l 


Hen 


14 




Cn 


54 


Eber& Butler (1988 1 


Cm 


3 




Ni 


3 


Becker & Butler {1989} 


Nn 


50 




Nm 


5 




Oi 


3 


Becker & Butler (19881 


On 


52 




Om 


3 




Mgi 


88 


Przybilla et al. EOoTl 


Mg ii 


37 




Aim 


12 


Dufton et al. fi986t 


Sin 


34 


Trundle et al. ( 2004 1, see also text 


Sin 


28 




Si iv 


18 




Sn 


78 


Vrancken et al. 119961 


Sm 


21 




Fen 


8 


Butler et al., in prep. 


Fein 


264 




Feiv 


16 





As will be shown in Sect. |4] the programme stars span a 
restricted range in T e ff and logg. This implies that the depar- 
tures from LTE should be roughly uniform among our sample 
and the relative abundances fairly insensitive to the treatment 
of NLTE line formation. Extremely few NLTE iron abundance 
studies of early B-type stars have been presented in the litera- 
ture (to our knowledge, the only one being the analysis of the 
B1.5 dwarf HD 35299 by Vrancken [TW71 . The magnitude of 
the NLTE corrections affecting the Fe m lines, Ae=(log 6)nlte- 
(log 6)lte, hence deserves further discussion. The departures 
from LTE are found to negligible for the coolest objects in our 
sample (Ae 5; +0.05 dex), but to significantly increase with r e fr 
(e.g. Ae=+0.21 dex for yS Cep and +0.33 dex for CMa, re- 
spectively). Previous LTE values in the literature derived from 
Fern lines for Bl dwarfs or earlier subtypes can thus be re- 
garded as strict lower limits. For the other elements, an illus- 
trative LTE abundance analysis carried out on (3 CMa suggests 
systematically negative mean NLTE corrections never exceed- 
ing 0.17 dex (see Table0}. Such modest departures from LTE 
for early B-type dwarfs or subgiants explains the reasonable 
agreement with past LTE abundance determinations in the lit- 
erature (e.g. Martin 2004 and Rolleston et al. 2003 in the case 
of y Peg and CMa, respectively); large discrepancies can 
generally be ascribed to other factors, such as differences in the 
adopted atmospheric parameters or atomic data. 



Table 4. NLTE corrections in the case of (3 CMa. The number 
of lines considered is given in brackets. Note that the listed 
corrections for some ions are based on very few lines and may 
not be representative. 



Element 


Ion 


Ae=(lo2 6)mitp— (I02 6)1 tp 


c 


C ii (8) 


-0.07 




Cm (1) 


0.00 




C n+C m (9) 


-0.06 


N 


Nn (25) 


-0.11 




Nni(l) 


0.00 




Nii+Nm(26) 


-0.11 


O 


On (30) 


-0.14 


Mg 


Mgn(2) 


-0.10 


Al 


Al m (4) 


0.00 


Si 


Sin (2) 


+0.03 




Sim (6) 


-0.26 




Siiv(l) 


-0.06 




Sin+Sim+Siiv (9) 


-0.17 


S 


Sn(l) 


+0.10 




Sm(l) 


-0.12 




Sn+Sm(2) 


-0.01 


Fe 


Fern (25) 


+0.09 



3.2. Estimate of atmospheric parameters 

The stellar physical parameters were derived in a standard way, 
whereby T e g and log g are self-consistently determined using 
an iterative scheme. Few iterations are needed before conver- 
gence is attained when sensible starting values can easily be 
estimated a priori, as is the case here. The effective tempera- 
ture was estimated from the Si ii/m/iv ionization equilibrium. 
Unfortunately, the spectral lines of these three ions are only 
measurable in ft CMa and ft Cep. For the other stars, we had to 
solely rely on the Si 11/111 or Si m/iv lines. This has motivated us 
to replace the 12-level Sin model ion initially implemented in 
DETAIL (Becker & Butler 1990) by a more complete atomic 
model with 34 levels developed by D. J. Lennon. This is ex- 
pected to allow for a more realistic treatment of the ioniza- 
tion balance and leads to a slightly cooler temperature scale 
(see Trundle & Lennon 2005 1. Figure [5] provides various ex- 
amples of calibrations between the Si line ratios and T e g. We 
used solar Si abundances (as appropriate for our stars; see be- 
low) and atmospheric models with a £ mo d value matching as 
closely as possible the microturbulence estimated from the line 
analysis. Several line ratios (up to 12 in total) were used for the 
temperature determination. Not all of them are completely in- 
dependent (only seven transitions are involved), but consistent 
results were obtained in all cases. The dispersion between the 
obtained values was taken as representative of the temperature 
uncertainty. This proved to be comparable from star to star and 
of the order of 1000 K. In cases where the abundances of C, N 
and S are based on lines of two adjacent ionization stages, no 
convincing evidence for systematic discrepancies was found. 
As a result, all abundances yielded by the individual lines were 
averaged irrespective of their ionization stage. This gives cre- 
dence to the reliability of our temperature scale, although the 
scope of such a comparison is limited by the coarse treatment 
of the minor ionic species C m and N m in our model atoms 
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Fig. 2. Examples of calibrations between various Si line ratios 
and the effective temperature, as a function of £. The depen- 
dence on the microturbulence must in certain cases be taken 
into account, as the lines used may be of very different strength. 
Solid line: £=5 km s _1 ; dotted line: £=10 km s ; dashed 
line: £=15 km s . In this case we adopted log g-3.5 [cgs], 
log e(Si)=7.20 dex and £ mo d=8 km s _1 . 



(Table QJ. We will comment further on this point in Sect. 14.11 
by comparing our T e ff values with literature data. 

The surface gravity was derived by fitting the collisionally- 
broadened wings of the Balmer lines. Line broadening mech- 
anisms such as rotation or microturbulence only affect the line 
cores and hence no prior knowledge of these quantities is 
needed. The Hfi line was not used in /3 Cep, as the Ha pro- 
file is clearly filled-in by a variable emission component. HyS 
may thus slightly be affected by emission as well. All of He, 
HS, Hy and H/3 were used for the other stars, as the lower se- 
ries members are unlikely to be substantially affected by cir- 
cumstellar/wind emission. Indeed, they systematically yielded 
consistent results, as can be seen in Fig.[5]in the case of ft CMa. 
Typical uncertainties are of the order of 0.15-0.20 dex and re- 
flect both the dispersion of the values obtained using the var- 
ious lines, errors arising from the uncertainty on the T e ff cali- 
bration and imperfections in the merging of the echelle orders. 
As mentioned in Sect. |3 this step in the reduction procedure 
was not possible for the McDonald data. Therefore, the profile 
fitting of the Balmer lines in V2052 Oph was performed using 
spectra obtained with the Coude Echelle spectrograph mounted 
on the 2-m telescope of the Tautenburg observatory, Germany 
(see Neiner et al., in prep.). For 12 Lac, we made use of a snap- 
shot spectrum extracted from the ELODIE archives. 

The microturbulence, £, was estimated by requiring the 
abundances yielded by the individual O n lines to be indepen- 
dent of the line strength. To estimate the uncertainties, we var- 
ied this parameter until the slope of the log e(Q)-\og(EW I A) 
relation differs from zero at the 3cr level. The oxygen features 
are the strongest metal lines in our spectra and are thus best 
suited for this purpose. However, the microturbulent veloci- 
ties obtained from the Nn lines are much lower than those 
obtained from the On features for the most massive objects 
(ft CMa, CMa, ft Cep and 12 Lac) and are in better agree- 
ment with the values inferred for B-type stars in general us- 
ing other methods (e.g. Fitzpatrick & Massa 2005 1. Although 
extensively discussed in the literature, the possible causes of 
these inconsistencies remain unclear (e.g. Trundle et al. 2004; 
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Fig. 3. Synthetic {dotted line; r eff =24000 K and logg= 
[cgs]) and observed (solid line) spectra of ft CMa for the 
gions encompassing the Balmer lines. 



3.5 
re- 



Table 5. Abundance differences when adopting the microtur- 
bulent velocities yielded by the N n lines (£=2 and 3 km s -1 for 
£' CMa and ft Cep, respectively) instead of the default values 
derived from the O n features (£=6 km s in both cases). 





f 1 CMa 


/3Cep 


A£ 


from 6 to 2 km s~' 


from 6 to 3 km s" 1 


AHe/H 


+0.014 


+0.005 


Aloge(C) 


+0.02 


+0.01 


A log e(N) 


+0.08 


+0.05 


Aloge(O) 


+0.20 


+0.13 


Aloge(Mg) 


+0.18 


+0.14 


A log e(Al) 


+0.10 


+0.06 


Aloge(Si) 


+0.09 


+0.17 


Aloge(S) 


+0.04 


+0.02 


A log e(Fe) 


+0.06 


+0.04 



see also Sect. 15. II . We provide in Table[5]some illustrative ex- 
amples of the abundance differences resulting from the exact 
choice of the species used to determine the microturbulence. 
Although we note that adopting the low values suggested by 
the N ii lines would lead to unreasonably high abundances for 
several strong lines of other elements, we cannot rule out a sys- 
tematic, albeit slight (A log e < 0.2 dex), underestimate of the 
stellar abundances in some objects. 

3.3. Estimate of chemical abundances 

The abundances are derived once the atmospheric parameters 
mentioned above are known by matching the measured EWs of 
the selected lines (tabulated in Table IA.fl and the values mea- 
sured in the synthetic spectra. Such an EW-based abundance 
analysis is made possible by the fact that all our stars are rel- 
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Table 6. Calculation of the error budget in the case of /3 CMa. 
o"i nt : line-to-line scatter; cr TcS : variation of the abundances for 
Ar eff =+1000 K; o- logg : as before, but for Alogg=+0.15 dex; 
erf as before, but for A£=+3 km s _1 ; ctt^i io Sg : as before, but 
for Ar e ff=+1000 K and Alogg=+0.15 dex; cr r : total uncer- 
tainty. 





0"tot 






°> 




<T 7 - 


AHe/H 


0.014 


0.010 


0.007 


0.020 


0.003 


0.028 


Aloge(C) 


0.037 


0.074 


0.001 


0.019 


0.059 


0.103 


Aloge(N) 


0.121 


0.046 


0.008 


0.026 


0.042 


0.139 


Aloge(O) 


0.088 


0.078 


0.040 


0.109 


0.053 


0.174 


Aloge(Mg) 


0.156 


0.060 


0.015 


0.050 


0.050 


0.182 


A log e(Al) 


0.096 


0.093 


0.022 


0.032 


0.060 


0.152 


Aloge(Si) 


0.206 


0.036 


0.009 


0.086 


0.028 


0.228 


Alog6(S) 


0.212 


0.100 


0.035 


0.015 


0.060 


0.245 


A log 6(Fe) 


0.178 


0.033 


0.031 


0.015 


0.046 


0.190 



atively slow rotators, as will be shown below. For consistency, 
direct integration was used in both cases. In a few instances, 
multicomponent fitting using a Voigt profile was performed 
with IRAF tasks in the case of well-separated blends. 

Five sources of error on the final abundances were consid- 
ered. Apart from the line-to-line scatter 4 (o"i nt ), we first cal- 
culated the errors arising from the uncertainties on the atmo- 
spheric parameters (cr rdT , o"i ogs and cr^). They were derived by 
computing the abundances using models with atmospheric pa- 
rameters deviating from the final values by the relevant uncer- 
tainties. We also took into account the fact that the determi- 
nations of r e ff and \ogg are strongly coupled (cr rdT /io g g)- The 
other covariance terms do not significantly contribute to the to- 
tal error budget and can be neglected. Finally, we quadratically 
add these errors to obtain the total uncertainty, cr T . Table [5] il- 
lustrates these calculations in the case of f3 CMa. It can be seen 
that the accuracy of the abundance determination varies signifi- 
cantly depending upon the element considered, being typically 
of the order of 0. 1 dex for C, but much larger for S and Si, for 
instance (~0.2-0.3 dex). This primarily reflects the sensitivity 
to changes in the physical parameters, limitations in the model 
atoms, as well as inaccuracies in the atomic data. 

A comment is necessary regarding the helium abundance 
determination in V2052 Oph. This star has been classified as 
He-strong by Neiner et al. (2003 1 who obtained He/H=0.21 
by number based on a detailed NLTE analysis of a number 
of optical He i lines. Although our curve-of-growth analysis 
does not support such a high helium content (see below), con- 
trary to the other stars in our sample we were unable to ob- 
tain consistent fits for all the He i lines considered. In partic- 
ular, it was impossible to fit the wings of the diffuse lines us- 
ing a single helium abundance, some lines being well fit with 
He/H=0.09, others with twice the solar value (particularly the 
singlets). Interestingly, a similar phenomenon is observed in f3 
Cep. A detailed investigation is beyond the scope of this paper, 
but this phenomenon may be related to the existence of surface 
inhomogeneities/vertical stratification of helium in these two 



4 The Mg abundance was often derived from a single line. In this 
case, we fixed cr mt to 0.15 dex. 



magnetic stars. As mentioned in Sect. 13. II Kurucz models with 
He/H=0.178 have been used for V2052 Oph. This choice has 
a negligible impact on the atmospheric parameters and abun- 
dances (AHe/H=0.003; A log e < 0.05 dex). 

3.4. Estimate of the amount of line broadening 

We provide in Table [7] a total line broadening parameter, vt, 
defined as: 

v\ = (OR sin if + v^ acro + v 2 puls , ( 1 ) 

with O the equatorial angular rotational frequency, R the stel- 
lar radius, i the inclination angle, v macro the macroturbulent ve- 
locity and Vp U i s the amount of pulsational line broadening av- 
eraged over the oscillation cycle. The latter is expected to be 
significant for non-radial modes, but to be small for radial pul- 
sators (cf. comparison between Tables[2and 0. The amount of 
macroturbulent broadening is unknown for our objects which 
are on, or close to, the main sequence. The vj values were 
estimated as a final step of the abundance analysis by fitting 
the profiles of some isolated O n lines with a grid of synthetic 
spectra with the appropriate microturbulent velocity and con- 
volved with a rotational broadening profile (Gray [T97rjl . as- 
suming a limb-darkening coefficient, e=0.4 (Claret 2000 1. The 
instrumental profile was estimated from calibration lamps. 

To disentangle the effect of rotational and macroturbulent 
broadening on the one hand, and of the averaged pulsational 
line broadening on the other hand, we have selected the spec- 
trum least affected by pulsations in the time series, i.e. with the 
narrowest and most symmetric line profiles (this is not possible 
for V836 Cen, as we only have one observation). These were 
used to derive the quantity: (OR sin i) 2 + v^ acro . The amount of 
averaged pulsational line broadening immediately follows from 
Eq. [0 The true stellar rotation rates, £1R, are known for most 
targets (Sect. l5~3l and indicate that some macroturbulence does 
seem to occur in some stars (see Table[7). 

4. Results and comparison with previous studies 

4.1. Physical parameters 

The atmospheric parameters of the programme stars are pre- 
sented in Table We also provide the revised MK spectral 
types based on our T e g estimates and the spectral type-r e ff cal- 
ibration for Galactic OB dwarfs of Crowther ( 1998 1. It is nec- 
essary to assess the reliability of our derived T e ff and log g esti- 
mates before discussing the chemical properties of our targets. 
An exhaustive literature search was not attempted. Instead, a 
comparison was made with atmospheric parameters obtained 
using three independent methods: (i) T e ff primarily derived 
from Stromgren photometry and logg from profile fitting of 
Hy (GL); (2) both parameters estimated from Walraven colour 
indices, or from Geneva and Stromgren photometry when these 
data were not available (Heynderickx et al. 1994); (iii) tempera- 
tures based on the continuum fitting of low -resolution IUE data 
(Niemczura & Daszyhska-Daszkiewicz 2005 1. As can be seen 
in Fig. 0] our temperatures (and consequently surface gravi- 
ties) are systematically lower than the values quoted by GL, but 
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Table 7. Physical parameters of the programme stars. The revised MK spectral types are based on our T e fi estimates and the 
calibration of Crowther ( 1998 ». The 1-cr uncertainty on r ei f and logg are 1000 K and 0.15 dex, respectively (except for V836 
Cen and V2052 Oph: A log g=0.20 dex). 





y Peg 


<5Cet 


vEri 


,8 CMa 


f 1 CMa 
Max EWs Min EWs 


V836 Cen 


V2052 Oph 


/3Cep 


12 Lac 


Revised spectral type 


B1.5-B2 V 


B1.5-B2 V 


B1.5-B2 V 


B1.5 V-IV 


B0.5-B 1 V-IV 


B1.5 V 


B1.5-B2 V 


B 1 Vevar 


B1.5 V 


T eB (K) 


22 500 


23 000 


23 500 


24000 


27 000 28 000 


24 500 


23 000 


26 000 


24 500 


iogr dr (K) 


4.352±0.019 


4.362±0.019 


4.371±0.019 


4.380±0.018 


4.431±0.016 4.447±0.016 


4.389±0.018 


4.362±0.019 


4.4 15 ±0.0 17 


4.389±0.018 


log g (cms 2 ) 


3.75 


3.80 


3.75 


3.50 


3.70 3.80 


3.95 


4.00 


3.70 


3.65 


f (kms- 1 ) 


l -i 




10±4 


14±3 


6±2 


6±3 


m 


6±3 


10±4 


vj- (km s _1 ) 


10 


14 


36 


23 


10 


15 


62 


29 


42 


[(nRsini) 2 + i-Lcro]" 2 (km s _1 ) 


-10 


-13 


-16 


-20 


-10 




-61 


-29 


-37 


Vpuk (km s -1 ) 


negligible 


negligible 


-32 


-11 


negligible 




negligible 


negligible 


-20 


OR (kms- 1 ) 




14 or 28° 


6" 


32* 




2" 


56' 


26 c 


45'" 


v macro (km S~ ] ) 




? 


> 15 


9 






> 20 


> 10 


? 



" From seismic studies (Dupret et al. 2004 Pamyatnykh et al. 2004 Aerts et al. 120061 . 

b From modelling of the line-profile variations (the so-called "moment method 1 ': Aerts 1996 Desmet et al. 200*61. 

c Calculated from the stellar radius derived from evolutionary tracks (Fig.|^J and assuming that the rotational period can be identified with the recurrence timescale of the changes 
affecting the UV lines: P rol =12.001 06 d in/3 Cep (Henrichs et al. l2TOI and 3.638 833 d in V2052 Oph (Neiner et al. l2*0*6*3*l . 



there are reasons to believe that their temperature scale is too 
hot. In particular, they applied a systematic upward correction 
to their temperatures in order to match the empirical calibration 
derived by Code et al. i 19761 1 from angular diameter and ab- 
solute flux measurements. However, more recent and accurate 
observations (Smalley & Dworetsky 1995 1 lead to systemati- 
cally lower values by about 1000 K in this temperature range 
(see discussion in Lyubimkov et al. 2002 1. The star /? CMa was 
incidentally included in this sample of standard stars and is ide- 
ally suited for a direct comparison: r e fr was decreased from 
25 180+1 1 30 (Co de et al. 119751 to 24020±1 150 K (Smalley & 
Dworetsky 11995b . which is now in better agreement with our 
results. Additionally, we note that Aufdenberg et al. ( 1999) ob- 
tained absolutely identical values to ours for r e fr and log g from 
the fitting of the spectral energy distribution from the near-UV 
to the infrared using a fully line-blanketed, spherical and NLTE 
atmospheric model. On the other hand, our T e s values tend to 
be slightly higher than those obtained from colour indices. The 
photometric calibrations in the Walraven photometric system 
are based on the grids of lightly line-blanketed model atmo- 
spheres of Kurucz dl979> . Unfortunately, the lack of further 
details (e.g. zero points) does not allow us to identify the cause 
of these discrepancies. Better agreement is found with the tem- 
peratures derived from UV data, but large differences are found 
for our two hottest objects: ^ CMa and /3 Cep. As regards the 
surface gravities, our values are identical, within the errors, to 
those presented by Heynderickx et al. ( 1994). 

A variation Ar e (f~1000 K and Alogg~0.1 dex is inferred 
during the radial pulsation cycle of CMa. This is com- 
parable to the uncertainties, but the changes in the EWs of 
temperature-sensitive lines (Table lA~Tt are such that the tem- 
perature variations are clearly real. Furthermore, Beeckmans 
& Burger ( 1977 1 obtained very similar results from UV colour 
indices: Ar er r= 1400+430 K. This can be compared to values 
obtained by various means for other /3 Cephei stars dominated 
by a large-amplitude radial mode: e.g., Ar e ff=4000±2000 K 
and Alogg=0.7±0.4 dex for cr Sco (van der Linden & Butler 
[T9g8l , Ar eff =900±120 K for V2052 Oph (Morto n & H ansen 
fT974l or A7W-1200 K for 6 Cet (De Ridder et al. l2002l 



4.2. Stellar abundances 

The mean NLTE abundances are given in Table along with 
the resulting metallicity, Z. As a consistency check, we have 
repeated the abundance analysis of CMa for the two phase 
ranges considered, finding nearly identical results in both cases. 
A detailed abundance analysis has been performed for seven 
stars in our sample by GL. However, a detailed comparison is 
probably not warranted for the following reasons: (a) the two 
sets of atmospheric parameters adopted are systematically dif- 
ferent, with their temperature scale being likely too hot (see dis- 
cussion above); (b) the NLTE corrections applied are based on 
atmospheric models with an inadequate treatment of metal line 
blanketing (Gold 1984). Nitrogen being of particular relevance 
(see below), we simply note the perhaps fortuitous close agree- 
ment between the NLTE abundances, with differences never 
exceeding 0.2 dex. Ignoring previous LTE studies (e.g. Watson 
119711 Peters 11973b . the same quantitative agreement is found 
for y Peg, 6 Cet and V2052 Oph with other NLTE CNO abun- 
dance determinations in the literature (Andrievsky et al. 1999, 
Korotin et al. ll999al b.c; Neiner et al. l2003t . 

5. Discussion 

5. 1 . Neglect of stellar winds 

It is natural to ask to what extent the results presented in this 
paper are robust against the neglect of the stellar wind, partic- 
ularly in view of the fact that the inferred microturbulent ve- 
locities are in some cases comparable to the sound velocity in 
the O ii line-forming regions (recall that an overestimate of this 
parameter would lead to spuriously low abundances). As ex- 
pected, relaxing the assumption of LTE leads to a decrease of 
this quantity, but relatively large values are still obtained with 
a full NLTE treatment. In the case of j3 CMa, for instance, £, 
decreases from 18 to only 14 km s -1 when NLTE is enforced. 
Furthermore, although our sample is admittedly small and the 
range in surface gravity limited, there is an indication for higher 
microturbulent velocities in the more evolved objects (see also 
Kilian ll992l and Daflon et al. l2004> . It is tempting to interpret 
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20 22 24 26 28 20 22 24 26 28 20 22 24 26 28 30 
T en ( GL ) LkKj T e(f (Heynderickx) LkKj T e(f (Niemczura) LkK] 

Fig. 4. Comparison between our T e g values and literature data. 

Table 8. Mean NLTE abundances, along with the total 1-cr uncertainties, crj (by convention, log e[H] = 12). The number of used 
lines is given in brackets. For comparison purposes, the last four columns give: the unweighted mean values for our sample 
(Mean), typical values found for OB dwarfs in the solar neighbourhood (Daflon & Cunha 2004 OB stars), the standard solar 
composition of Grevesse & Sauval ( 1998 Sun 1-D), and finally solar abundances recently derived from 3-D hydrodynamical 
models (Asplund et al. 2005 Sun 3-D). The metallicity, Z, is calculated assuming that the abundances of the elements not under 
study is solar (Grevesse & Sauval fT9"9"%l . We define [N/C] and [N/O] as log[e(N)/e(Q] and log[e(N)/e(0)], respectively. 





7 Peg 


<5Cet 


v Eri 


/3CMa 


f 1 


CMa 


V836 Cen 


V2052 Oph 












Max EWs 


Min EWs 






He/H 


0.079±0.()25 (9) 


0.096±0.024(10) 


0.076±0.024 (8) 


0.070±0.028 (9) 


0.098±0.017 (10) 


0.099*0.016 (10) 


0.078±0.024(10) 


0.118±0.032 (9) 


loge(C) 


8.20±0.05 (8) 


8.09±0.08 (6) 


8.24±0.12(10) 


8.16±0.11 (9) 


8.19±0.11 (9) 


8.17±0.12(9) 


8.34±0.13 (8) 


8.21±0.07 (4) 


loge(N) 


7.58±0.11 (23) 


8.05±0.11 (26) 


7.87±0.09 (18) 


7.59±0.14(26) 


8.00±0.16(34) 


7.99±0.18(34) 


7.73±0.10(21) 


7.99±0.17 (10) 


logf(O) 


8.43±0.28 (21) 


8.45±0.26 (19) 


8.51±0.24(19) 


8.62±0.18(30) 


8.59±0.16(34) 


8.59±0.17 (34) 


8.49±0.24 (25) 


8.39±0.30(14) 


log £(Mg) 


7.44±0.25 (1) 


7.52±0.29 (1) 


7.38±0.24(1) 


7.30±0.18(2) 


7.37±0.20(1) 


7.38±0.20(1) 


7.41±0.19(3) 


7.35±0.32(1) 


loge(Al) 


6.12±0.16(4) 


6.13±0.21 (4) 


6.08±0.16(4) 


6.00±0.15(4) 


6.16±0.21(4) 


6.16±0.23(4) 


6.13±0.11 (4) 


6.09±0.24 (2) 


loge(Si) 


7.19±0.29(6) 


7.28±0.29 (6) 


7.21±0.26 (9) 


7.17±0.23 (9) 


7.14±0.23(4) 


7.12±0.19(4) 


7.14±0.19(9) 


7.16±0.37 (6) 


logf(S) 


7.22±0.20 (14) 


7.26±0.24 (4) 


7.32±0.22 (5) 


7.14±0.25 (2) 


6.98±0.11 (2) 


7.00±0.21 (2) 


7.39±0.26 (2) 


7.32±0.21 (5) 


log f(Fe) 


7.25±0.16(28) 


7.32±0.18(21) 


7.36±0.19(17) 


7.17±0.19(25) 


7.31 ±0.19 (32) 


7.28±0.25 (32) 


7.30±0.09(16) 


7.37±0.21 (9) 


Z 


0.0091 ±0.0021 


0.0100±0.0021 


0.0105±0.0022 


0.0104±0.0021 


0.0110±0.0018 


0.0108±0.0020 


0.0105±0.0022 


0.0097±0.0022 


[N/C] 


-0.62±0.12 


-0.04±0.14 


-0.37±0.15 


-0.57±0.18 


-0.19±0.2() 


-0.18±0.22 


-0.61±0.17 


-0.22±0.19 


[N/O] 


-0.85±0.30 


-0.40±0.29 


-0.64±0.26 


-1.03±0.23 


-0.59±0.23 


-0.60±0.25 


-0.76±0.26 


-0.40±0.35 





/JCep 


12 Lac 


Mean 


OB stars 


Sun 1-D 


Sun 3-D 


He/H 


0.078±0.028 (8) 


0.075±0.()31 (8) 


0.085±0.015 


0.10" 


0.085 


0.085 


log 6(C) 


8.02±0.10(11) 


8.22±0.12(6) 


8.18±0.09 


-8.2 


8.52 


8.39 


loge(N) 


7.91±0.13 (19) 


7.64±0.18(17) 


7.82±0.19 


-7.6 


7.92 


7.78 


log 6(0) 


8.47±0.14(30) 


8.42±0.23 (22) 


8.49±0.08 


-8.5 


8.83 


8.66 


loge(Mg) 


7.31±0.21 (1) 


7.29±0.23 (1) 


7.37±0.07 


-7.4 


7.58 


7.53 


loge(Al) 


6.02±0.16(4) 


6. 11 ±0.17 (3) 


6.09±0.05 


-6.1 


6.47 


6.37 


loge(Si) 


7.11±0.23 (8) 


7.11±0.27(6) 


7.17±0.05 


-7.2 


7.55 


7.51 


loge(S) 


7.14±0.37 (2) 


7.10±0.31 (2) 


7.21±0.13 


-7.2 


7.33 


7.14 


log e(Fe) 


7.24±0.23 (23) 


7.30±0.20 (22) 


7.29±0.06 


-7.4* 


7.50 


7.45 


Z 


0.0091 ±0.00 13 


().0089±0.0018 


().0099±0.0007 


-0.0099 


0.0172 


0.0124 


[N/C] 


-0.11±0.17 


-0.58±0.22 


-0.36±0.21 


—0.6 


-0.60 


-0.61 


[N/O] 


-0.56±0.19 


-0.78±0.29 


-0.67±0.21 


—0.9 


-0.91 


-0.88 



a From Lyubimkov et al. 120041 . The primordial helium abundance is He/H-0.08 (e.g. Olive & Skillman 2004 . 
h Mean NLTE abundance of four BA supergiants in the solar vicinity (Przybilla et al. 120061 . 



these high microturbulences as an artefact of our assumption 
of a hydrostatic photosphere (e.g. Lamers & Achmad 1994 1. 
However, there is ample evidence in the literature indicating 
that high microturbulent velocities are needed, even when state- 
of-the-art, "unified" codes accounting for mass loss and spher- 
ical extension are used (e.g. Mokiem et al. 2005 1. 

A number of studies have investigated in recent years the 
differences between the atmospheric parameters and/or abun- 
dances of early-type stars obtained using the plane-parallel 
code TLUSTY (Hubeny & Lanz ll99"5l and the unified codes 
CMFG EN/FASTWIND (Hillier & Miller Il998l Puis et al. 
2005 1. This comparison was performed for a small sample 
of Galactic/SMC B-type supergiants (Urbaneja et al. 2003 



Trundle et al. 2004 Dufton et al. 120051 and Galactic O-type 
dwarfs (Bouret et al. 2005 , Martins et al. 2005 1. The differences 
were claimed to be of the order of the uncertainties and hence 
of no significance. These consistency checks support the valid- 
ity of our analysis, especially considering that weaker winds 
by at least two orders of magnitude are expected in our sample. 
Perhaps of more relevance is the fact that the same conclusion 
was reached for some late O-type dwarfs in the SMC having 
mass-loss rates comparable to what is anticipated for/3 Cephei 
stars (M down to 10 9 -10 10 M Q yr _1 ; Bouret et al. l2003l . A 
very weak wind was derived for the most evolved star in our 
sample, ft CMa, from soft X-ray observations: M-6 x 10 9 
M yr _1 (Drew et al. [l994l 
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In summary, although forthcoming analyses of our pro- 
gramme stars using unified models would obviously be desir- 
able, the choice of the line formation code used is very unlikely 
to affect our conclusions about the abundance patterns and does 
not alleviate our concerns regarding the high microturbulent 
velocities. 

5.2. Abundance patterns 

As can be seen in Table [8] the abundances of He, C, O, Mg, 
Al, Si and Fe span a limited range in the studied stars (less 
than 0.3 dex). This is comparable to the total uncertainties and 
indicates that the abundances of these species are remarkably 
uniform among our sample. The homogeneous nature of our 
sample in terms of spectral type implies that our results should 
be fairly insensitive to several sources of systematic errors. 
Nevertheless, we verified that no significant trends exist be- 
tween the abundances and T e ff . This important check supports 
the reliability of our NLTE corrections and temperature scale 
and in turn partly explains the small scatter within our sam- 
ple. Nitrogen and sulphur exhibit a significantly larger spread 
of the order of 0.4-0.5 dex. The S abundances are derived from 
very few lines (sometimes of different ionization stages), which 
suggests that this scatter may not have a physical basis. On the 
contrary, we will argue below that the star-to-star N variations 
are real and may be interpreted in terms of mixing processes 
(Sect.|5J}- 

One issue we wish to address with the present study is the 
possible existence of a dichotomy between the chemical prop- 
erties of the f3 Cephei stars and those of non-pulsating early 
B-type stars. In the following we shall use the mean abundance 
of OB stars in the solar neighbourhood as a baseline for com- 
parison purposes. The final results of a comprehensive, homo- 
geneous abundance study of about 90 late O- to early B-type 
Galactic stars on (or close to) the main sequence have recently 
been presented by Daflon & Cunha ( 2004 1. The analysis is sim- 
ilar in several aspects to our study (use of DETAIL/SURFACE 
and Kurucz models), a fact which reduces the systematic errors. 
The typical abundances at the solar galactocentric distance (all 
our stars lie within 1 kpc) are quoted for all the elements in 
Tabled along with the mean values for our sample. 5 A direct 
comparison shows a near-perfect agreement for He, C, O, Mg, 
Al, Si and S. Our mean LTE iron abundance (~7. 18 dex) is in 
reasonable agreement with the value of ~7.3 dex estimated by 
Daflon et al. (1 19991 1200 Tat for early B-type dwarfs. However, 
we prefer to compare our results to the NLTE Fe abundances of 
four relatively nearby (d % 3 kpc) late-B to early-A supergiants 
obtained by Przybilla et al. (2006 1 using a set of optical Fen 
lines and the model atom developed by Becker ( 1998 1. The Fe 
abundance is not altered during the evolution of massive stars, 
allowing for a meaningful comparison. Our mean NLTE abun- 
dance is slightly lower than the value quoted by Przybilla et 



5 It is customary in seismic studies to use the standard solar mix- 
tures of Grevesse & Sauval 1 1998 1 or more recent values based on 3- 
D hydrodynamical models, as summarized by Asplund et al. 1 2005 1. 
These values are also quoted for convenience in Table|S| 



al. (|2006), but is still compatible when taking the errors into 
account (Table [8}. 



5.3. Evidence for deep mixing 

As stated above, nitrogen is the chemical species exhibiting the 
largest star-to-star spread. This may indicate the existence in 
some stars of the products of (incomplete) CNO-cycle burning 
brought to the surface by deep internal mixing. Evolutionary 
effects are known to greatly enhance the photospheric N abun- 
dance in BA supergiants, but fast rotation leading to meridional 
currents and shear instabilities in the interior is widely believed 
to be responsible for mixing prior to the first dredge-up in sin- 
gle early-type stars. The positions of the programme stars in the 
HR diagram are displayed in Fig. [5] This shows that our targets 
are largely unevolved and have masses ranging from about 9 to 
20 M0. The most recent evolutionary models including rota- 
tion suggest that fast-rotating B stars in this mass range could 
display a measurable N excess of the order of a few tenths of a 
dex at the end of the main sequence. This should be accompa- 
nied by a slight C depletion at the 0. 1 dex level, which is com- 
parable to our uncertainties and therefore barely detectable. On 
the other hand, He and O should remain virtually unaffected 
(Heger & Langer l2000l Meynet & Maeder l2003l . The signa- 
ture of deep mixing in the stellar interior can be best revealed 
by examining the abundance ratios between N, C and O ([N/C] 
and [N/O]; see TableEJl. This shows that y Peg, /? CMa, V836 
Cen and 12 Lac have ratios typical of both nearby OB stars 
and the Sun, whereas 6 Cet, CMa, V2052 Oph, /3 Cep and 
to a much lesser extent v Eri display higher values by up to 0.6 
dex (a factor 4). The two subsamples of N-normal and N-rich 
stars are separated by about 0.5 dex in terms of the [N/C] ra- 
tio (which is the most robust diagnostic for an N enrichment), 
i.e. about three times the typical uncertainty in the determina- 
tion of this quantity (Tabled. The exact choice of the chemical 
species used to derive the microturbulence (oxygen or nitrogen) 
has little impact on these conclusions, as the [N/C] and [N/O] 
ratios are affected by at most 0.12 dex in the case of CMa 
and f3 Cep (see Table[5j. 

Such abundance patterns are in line with the predictions 
of the theoretical models cited above and support the idea that 
some objects are indeed experiencing substantial rotationally- 
induced mixing. Boron abundances are available for seven stars 
in our sample and can be used to further examine the relevance 
of this interpretation (ProfRtt & Ouiglev l2"0"0"Tl Venn et al.|2002 ; 
Mendel et al. 2006). This fragile element is destroyed by pro- 
ton capture at relatively low temperatures (T ~ 4 x 10 6 K) and 
is already completely depleted throughout most of the stellar 
envelope during the very early phases of evolution off the zero- 
age main sequence (ZAMS). The boron surviving in the outer 
stellar layers will quickly be transported downwards and be de- 
stroyed at the onset of rotational mixing. This will be followed 
by an increase of the nitrogen content as core -processed gas 
is being gradually mixed to the surface. As a result, one ex- 
pects the N-rich stars to be most severely B depleted and, con- 
versely, the stars with little B depletion to have near-solar N 
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Fig. 6. [N/C] and [N/O] ratios as a function of the NLTE B abundances. The boron data are primarily taken from Proffitt & 
Quigley (2001 1, but have been complemented in the case of 5 Cet, (3 Cep and 12 Lac by more stringent estimates from the 
literature (Venn et al. 2002 Mendel et al. 2006 1. The arrows denote upper limits. The abundance data are compared to the 



theoretical predictions of Heger & Langer (2000 1 for a 12 Mq star and three different values of the equatorial rotational velocity 
on the ZAMS: £1Rzams=99 (left-hand panels), 206 {middle panels) and 328 km s _1 {right-hand panels). The locus in each panel 
{dashed line and open squares) defines an age sequence with the time elapsed from the ZAMS increasing leftwards: from /=0 
to 15 Myrs (0 to 12.5 Myrs for £2fi ZA MS=99 km s~') by steps of 2.5 Myrs (see upper-right panel). The initial [N/C], [N/O], and 
boron abundances at f=0 have been taken as: -0.6, -0.9 (Table|S} and log e(B)=2.78+0.04 (meteoritic value; Zhai & Shaw 1994 1, 
respectively. 



abundances. This is consistent with the situation depicted in 
Fig. El 

However, high initial rotational velocities are required to 
reproduce the observed abundance patterns. For instance, an 
initial rotational velocity on the ZAMS, £1Rzams, of about 200 
km s is necessary to account for the N excess in 5 Cet (with an 
age of about 14 Myrs according to evolutionary tracks; Fig. [5}- 
The loss of angular momentum along the main sequence leads 
to a decrease of the rotation rate, but Q.R is still expected to 
be in this case of the order of 150 km s at the end of core 
H burning (Heger & Langer 120001 Meynet & Maeder l2003l . 
This is in excellent agreement with the observations of early-B 
dwarfs in young open clusters, which reveal only a slow decline 
of the rotation rate from -170 km s _I on the ZAMS to -130 
km s _1 at the end of core-hydrogen burning (Huang & Gies 
2006 1. Most stars enriched in nitrogen have just evolved off 
the main sequence (Fig.|5J, such that reproducing the observed 
abundance ratios requires high rotational velocities. This is in 
stark contrast with the fact that none of our targets is apparently 
rapidly rotating {Q.R sin i < 65 km s _1 ; TableQ. 

Comparisons between the surface abundances of the prod- 
ucts of nucleosynthesis and the predictions of models account- 
ing for mass loss and rotation are usually performed in a statis- 
tical sense (e.g. Herrero & Lennon 2004) and rely on the deter- 
mination of CNO abundances for a large sample of stars only 
having an estimate of their projected rotational velocity. One 
of the key features of this study is that the true rotation rates, 
OR, have been estimated for several of our targets by a detailed 



modelling of either the line-profile variations or the pulsation 
spectrum (see Tabled), offering an opportunity to constrain the 
models in a more quantitative way. First, seismic modelling of 
v Eri and V836 Cen has revealed a non-rigid rotation profile 
with very low equatorial velocities of about 6 and 2 km s _1 , re- 
spectively (Dupret et al. 2004 Pamyatnykh et al. l2004> . MOST 
observations of 6 Cet also lead to a likely rotational velocity of 
about 28 km s~' , or half this value depending on the exact mode 
identification (Aerts et al. 2006 1. Second, modelling of the line- 
profile variations yields QR-32 and 45 km s~' in [5 CMa and 
12 Lac, respectively (Desmet et al. 2006, Aerts 1996 1. In addi- 
tion, assuming that the periodicity of the changes affecting the 
UV lines can be identified with the rotational period (Neiner 
et al. l2003l Henrichs et al.|2000t, leads to QR-56 and 26 km 
s~' in V2052 Oph and /? Cep, respectively. Although this quan- 
tity is unfortunately unknown for CMa, this shows that at 
least three N-rich stars {6 Cet, V2052 Oph and /3 Cep) are in- 
trinsically slow rotators and should not exhibit any abundance 
peculiarities according to evolutionary models. 

Meridional currents are believed to be the most important 
mechanism governing the loss of angular momentum in (non- 
magnetic) B stars experiencing little mass loss, while shear 
mixing arising from differential rotation largely controls the 
transport of the chemical elements (Meynet & Maeder 2000 1. 
The observed discrepancy between the observations and the 
theoretical predictions could imply that the treatment of these 
two physical processes is still approximate in evolutionary 
models including rotation. It is also conceivable, for instance, 
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Fig. 5. Position of the programme stars in the HR diagram: (a) 
CMa, (b) B CMa, (c) B Cep, (d) 12 Lac, (e) v Eri, if) 5 Cet, 
(g) j Peg, (h) V2052 Oph and (0 V836 Cen. The luminosities 
were computed using Hipparcos parallaxes and the bolomet- 
ric corrections of Flower ( 1996 1. The extinction in the V band, 
Ay, was derived from the theoretical (B - V) colour indices of 
Bessell et al. (1998 1. The evolutionary tracks for solar metal- 
licity and without rotation are taken from Schaller et al. (1992 1. 
The initial masses are indicated on the right-hand side of this 
figure. The ZAMS is shown as a dashed line. 

that the internal rotation law, Q(r), is steeper than currently 
assumed. A decline of the rotational velocity by a factor 3-5 
from the core to the surface has been established from astero- 
seismological studies in two of our targets: v Eri (Pamyatnykh 
et al. l2004l and V836 Cen (Dupret et al. 2004). This appears 
broadly consistent with the predictions of theoretical models, 
but this piece of information is unfortunately lacking for the N- 
rich stars. Additional uncertainties arise from our very limited 
knowledge of the redistribution of angular momentum in mas- 
sive pulsating stars (see, e.g. Ando 1983, Lee & Saio 119931 . 
From the observational side, there appears to be an excess of 
slow rotators among the B Cephei class compared to the global 
B dwarf population, but an interpretation of this result is made 
difficult by the fact that stars with low OR sin i values tend to 
exhibit photometric changes of larger amplitudes and are hence 
more easily recognizable as B Cephei variables (Stankov & 
HandlerEHS}. 

It may be regarded as suggestive that 3 out of the 4 stars 
with an N excess possess a magnetic field with a longitudinal 
strength of -300 G C Ma; Hu brig et al. l2006l and -100 G 
(V2052 Oph: Neiner et al. l2003l B Cep: Henrichs et al. l2000l 
whereas similar observations for 3 stars with normal abun- 
dances (v Eri, B CMa and V836 Cen) did not yield any detec- 
tion (Schnerr et al. 120061 Hubrig et al. 120061 . No field has been 
so far reported in 6 Cet (Rudy & Kemp|^78 ), but much more 
sensitive observations are being planned (Briquet, private com- 



munication). In the models of Maeder & Meynet ( 2005 ), mag- 
netic fields tend to enforce near-solid body rotation and there- 
fore to suppress shear mixing. However, meridional circulation 
is in this case even more efficient at transporting chemical el- 
ements from the interior to the surface. Globally, one expects 
higher N surface abundances, but also the star to rotate faster. 
The inclusion of magnetic fields does not seem at first glance 
able to account for the observational data, but the models are 
still in their infancy and the theoretical predictions sensitive 
to some ill-defined quantities (e.g. diffusion coefficients). In 
this regard, we note that the predicted surface enrichments may 
strongly differ depending on the model considered (Maeder & 
Meynet 2005 ; Heger et al. 2005 1. Furthermore, dynamo action 
in the radiative layers is assumed to be driven by the so-called 
Tayler-Spruit instabilities (Spruit 2002 1, the existence of which 
is not firmly established and, in any case, poorly constrained 
from an observational point of view. Perhaps more importantly, 
complex field geometries are expected in the framework of this 
model, whereas V2052 Oph and B Cep possess a global dipole 
field (Neiner et al. 2003, Henrichs et al. 2000). This seems to 
argue for a magnetic field of fossil origin at least in these two 
objects. While it might be possible to explain in this context the 
existence of slowly-rotating magnetic B-type dwarfs depending 
on their pre-main sequence history (St§ pieh 2002 ), a nitrogen 
enrichment at the surface of such objects remains unexpected. 
Alternatively, one may postulate that these stars have entered 
the ZAMS with high rotational velocities, but have experienced 
a dramatic loss of angular momentum along the main sequence 
because of magnetic braking. However, this interpretation is 
very unlikely (see Donati et al. 2001 in the case of B Cep). 

One issue which needs to be clarified is whether the ob- 
served abundance ratios are only intrinsic to the B Cephei stars 
or are a characteristic feature of early B-type dwarfs in general. 
We show in Fig. the NLTE [N/C] and [N/O] ratios for our 
stars and 34 nearby B1-B2 dwarfs, as a function of Q/?sinz. 
When gathering the literature data, we restricted ourselves to 
stars spanning the same T e ff and log g ranges as our programme 
stars (see Table Q. In addition, only stars lying at less than 
about 2 kpc from the Sun were selected to avoid possible [N/C] 
and [N/O] variations arising from the chemical evolution of the 
Galaxy (Rolleston et al. 2000). It turns out that as many as 10 
stars in this sample are confirmed B Cephei stars (Stankov & 
Handler 2005 ). They are plotted, along with our targets, with 
different symbols in Fig.0 Great care must be exercised when 
interpreting these data in view of their heterogeneity, although 
using line ratios minimizes the systematic errors. The magni- 
tude of this problem can be assessed by comparing the [N/C] 
and [N/O] ratios for the two stars with independent determi- 
nations: HP 61068 (GL; Kilian fT992l and HD 216916 (GL; 
Daflon et al. 2001ai. The difference is only 0.08 dex on aver- 
age, but can reach up to 0.19 dex. In spite of the limitations 
mentioned in Sect. 14.11 we can also compare our results with 
GL for the seven stars in common. Once again, the differences 
are small on average (0.12 dex), but are of the order of 0.3 
dex in the most extreme case. Table [8] suggests significantly 
higher mean [N/C] and [N/O] ratios by about 0.25 dex for our 
subsample of B Cephei stars compared to nearby, presumably 
non-pulsating B dwarfs (Daflon & Cunha 2004). However, this 
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is only marginally the case when considering all the data col- 
lected from the various sources in the literature shown in Fig. [7] 
(~0.1 dex). This illustrates the fact that comparing the chem- 
ical properties of these two populations at these levels is sen- 
sitive to the exact definition of the samples (e.g. evolutionary 
status) and, most importantly, to systematics in the abundance 
determinations. Confidently establishing a possible dichotomy 
between the abundance patterns must await the homogeneous 
analysis of a large sample. Another more subtle difficulty lies 
in the fact that stars classified as non pulsating might prove 
to be /3 Cephei stars when observed under closer scrutiny (see 
Telting et al. 2006). We thus conclude at this stage that there 
is no convincing evidence for a higher amount of mixing in p 
Cephei stars compared to the global B dwarf population. 

We note that only the radial pulsator CMa or the stars 
dominated by a radial mode, 6 Cet, V2052 Oph and p Cep, are 
N-rich. To further examine the potential role played by pulsa- 
tions, we have investigated a possible link between the N ex- 
cesses and several pulsational observables (e.g. velocity am- 
plitudes, oscillation frequencies), but found no correlation be- 
tween these quantities. 

6. Conclusion and prospects 

We have presented a detailed, self-consistent determination 
of the atmospheric parameters and NLTE abundances of nine 
prototypical p Cephei stars based on very-high quality spec- 
troscopic data covering the entire oscillation cycle in most 
cases. This study reveals that the abundance patterns are re- 
markably similar among our sample, with the exception of 
nitrogen which is enhanced by a factor 2—4 in four targets. 
Caution should therefore be exercised when considering appar- 
ently slowly-rotating, early B-type dwarfs as good tracers of 
the present-day nitrogen abundance of the interstellar medium, 
either in the Galaxy or in the Magellanic Clouds. We inter- 
pret this nitrogen enrichment, which is systematically accom- 
panied by a boron depletion, as the result of internal mixing. 
The existence of stars with normal nitrogen abundances, but 
with a strong boron depletion (e.g. 12 Lac; Fig. |6j) supports 
this conclusion (Fliegner et al. 1996). Mass transfer processes 
in close binaries can mimic such abundance peculiarities (e.g. 
Vanbeveren 1989) and/or to lead to a dramatic spin down be- 
cause of tidal effects (Huang & Gies 2006 ). However, only 
two stars in our sample are confirmed members of a (wide) 
binary system: y Peg ( l P ol b~370 yr; Chapellier et al. 2006) and 
P Cep (Porb~90 yr; Pigulski & Boratyn[T332}. As y Peg is not 
N-enriched, this argues against a binary origin. Furthermore, 
the other stars have all been the subject of long-term, exten- 
sive spectroscopic campaigns (except V836 Cen) and do not 
present any signs of binarity. For the sake of completeness, 
we have also computed the peculiar space velocities, v pec , of 
our targets from Hipparcos data and found that none of the N- 
enriched stars exceeds the canonical threshold for a runaway 
status (v pec t 40 km s _1 ; Blaauw 1961 ). Hence there is no evi- 
dence for a dynamical kick imparted by the previous explosion 
of a putative companion as a supernova. Alternatively, the stars 
could have formed out of interstellar material contaminated by 
the ejecta of Type II supernovae. This hypothesis can unfor- 



tunately not be tested by comparing the derived abundances 
with those of OB stars in their immediate vicinity, as none of 
the four stars is known to be part of an open cluster (Stankov 
& Handler 2005). Likewise, a search for nearby open clusters 
with abundance data for early-type members in the WEBDA 6 
database proved inconclusive. Nevertheless, such a scenario is 
unable to explain the boron data and seems very implausible. 

At least three N-rich stars {6 Cet, p Cep and V2052 Oph) 
are intrinsically slow rotators and should not expose such a 
large amount of core-processed material at their surfaces ac- 
cording to the most recent evolutionary models accounting for 
rotation. It should be emphasized that the reliability of the the- 
oretical models has so far mainly been tested from observations 
of OBA supergiants in the Magellanic Clouds (e.g. Trundle 
& Lennon 2005 Heap et al. 2006 for some recent examples). 
Studies of B dwarfs, such as the present work, are free of sev- 
eral uncertainties which may render a confrontation between 
theory and observations difficult for such evolved objects: (a) 
the uncertainties on the mass-loss rates and on the related loss 
of angular momentum as the star evolves, (b) their ambigu- 
ous evolutionary status (pre- or post-red supergiant phase), and 
(c) the difficulties in estimating the rotation rate of supergiants 
(e.g. Dufton et al. 2006 ). One way of reconciling the model pre- 
dictions with the observations could be to assume that the £2(r) 
law in the interior is much steeper than hitherto assumed. The 
existence of a magnetic field in most N-rich stars (Sect. 15. 3> 
indicates that theoretical models including the effect of mag- 
netic fields are probably more relevant for interpreting our ob- 
servations. A direct comparison with models assuming that the 
field is produced via dynamo action in the radiative layers (e.g. 
Maeder & Meynet 2005 ) may not be appropriate, however, as 
the fields might be of fossil origin in these objects. 

We have determined the abundances of all the chemical 
species contributing significantly to the global metallicity, Z 
(except Ne). It is of importance to emphasize that the mild N 
enrichment observed in some stars has little incidence on the 
metal content, as this element is only a minor contributor to Z 
(~5%). As can be seen in Tabled the metallicity of the studied 
P Cephei stars is unremarkable and typical of the values ob- 
served for early B-type dwarfs in the solar neighbourhood, in 
accordance with recent works based on UV data (Niemczura & 
Daszyhska-Daszkiewicz 2005|). It can be noticed that the Sun 
is markedly more metal-rich than young, nearby B stars. This 
long-standing problem still lacks a clear explanation and may 
constitute a major source of uncertainty when choosing the ap- 
propriate metallicity value in oscillation codes. 

Recent multisite, ground-based campaigns have brought a 
wealth of new high-quality data on p Cephei stars (e.g. Aerts 
et al. I2004al Handler et al. I2006i . and much more is ex- 
pected in the near future thanks to a plethora of space mis- 
sions (e.g. MOST, COROT, Kepler). The diagnostic power of 
seismic modelling of early-type stars will hence become enor- 
mous, but evolution and oscillation codes unfortunately have 
difficulties in keeping pace with these dramatic observational 
developments. In particular, concerns have recently been raised 

6 The Database for Galactic Open Clusters (WEBDA) is available 
online at: http : / /www.univie . ac . at/webda/ . 
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Fig. 7. [N/C] and [N/O] ratios as a function of Q/?sin/ for our programme stars (circles) and a sample of 34 nearby B1-B2 
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( 1998 1, diamonds: Kilian ( 1992) and Kilian et al. ( 1994), hexagons: Mathys et al. (2002). The filled geometrical symbols denote 
known ft Cephei stars (Stankov & Handler 2005 1. The dashed lines show the typical ratios in the absence of mixing (Table[8). 
The Q.R sin i values from Table0have been used for our programme stars. The colour coding in the online version of this journal 
is the following: red: programme stars, blue: literature data for other ft Cephei stars, black: presumably non-pulsating B stars. 



regarding the inability of standard models to reproduce the rich 
oscillation spectrum of 12 Lac and v Eri. In both cases, the in- 
corporation of gravitational settling and radiative levitation has 
been claimed as an ad hoc solution to help resolve this problem 
(Pamyatnykh et al. 2004; Handler et al. 2006J. We find that both 
the [N/C] and the [N/O] ratios are about 0.25 dex above solar 
in y Eri (Table|S|l, a star whose rotation rate only decreases by 
a factor three between the core and the surface (Pamyatnykh et 
al. 2004). This is likely significant and provides some indica- 
tion for deep mixing. Such an interpretation would challenge 
the existence of diffusion effects at least in this star, as such 
microscopic processes are known to be severely hindered, or 
even completely inhibited, by large-scale fluid motions in the 
interior. However, P.-O. Bourge and collaborators have called 
our attention to the fact that radiatively-driven microscopic dif- 
fusion, and not rotational mixing, could actually be responsible 
for the N-excess in ft Cephei stars. They have already shown 
that this process can indeed lead to the accumulation of iron in 
the transition region of ft Cephei stars, thus explaining the ex- 
citation of all the observed modes in 12 Lac and v Eri (Bourge 
& Alecian 2006 Bourge et al. 2006a). Their most recent re- 
sults show that radiative forces on nitrogen are always greater 
than those on carbon and oxygen near the surface, thus leading 
to separation and variable excess of N depending on the stel- 
lar parameters (Bourge et al. 2006b ). This intriguing possibility 
opens a new avenue for the interpretation of the abundance pe- 
culiarities discussed in this paper. 
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Appendix A: EW measurements 

Table |A~T1 lists the measured EWs for the studied stars. Our 
EWs are systematically larger than the values quoted by GL 
for the seven stars in common, but a similar trend is evident 
when comparing their values with previous measurements in 
the literature (see their Fig. 4). This is likely to result from their 
use of Gaussian fitting instead of the direct integration used in 
our case. The same conclusion holds for the EWs of £ CMa 
presented by Hambly et al. JT996). 
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Table A.l. EW measurements in mA. A blank may have various causes: the EW was not reliably measurable (spectral features 
lying in the wings of the Balmer and diffuse He i lines were avoided); the line fell into echelle inter-order gaps; was considered 
blended for the relevant temperature range; or yielded a discrepant abundance. The latter point notably applies to the C n /14267 
line which is notoriously known to yield spuriously low abundances (e.g. Lennon et al. 2003 although Nieva & Przvbilla l2006l 
appear to have recently solved this problem). 
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20.9 
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Max EWs 
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V2052 Oph 
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O ii (Continued) 






















5190.56 


7.2 








28.5 


26.7 










5206.71 








57.8 


51.8 


49.7 


30.9 




42.8 




6640.99 






33.7 


62.1 


61.3 


57.2 






47.7 


36.8 
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49.2 


95.9 


84.7 


78.0 
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4139.35 


17.6 
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5.3 


6.0 


















4164.82 








66.2 


56. 1 


51.1 


58.4 




51.7 


69.0 


4166.84 


13.6 








16.1 


14.3 


20.7 




18.5 


22.2 


4222.27 




21.5 


31.1 
















4248.77 


4.8 




7.7 
















4261.39 


5.3 


5.1 






4.5 


4.3 




3.6 


4.2 


5.7 


4296.85 


11.3 






23.6 


19.9 


16.8 




12.1 


20.0 


21.9 


4310.35 


15.2 




36.7 


34.7 


27.8 


23.3 


27.6 


14.0 


21.4 


27.4 


4419.60 


32.7 


34.2 


52.4 


44.0 


30.6 


22.5 


44.8 


27.8 
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5.2 
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11.4 
















5833.94 






79.2 


80.1 
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11.9 










12.1 
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6.5 


9.2 






11.6 


9.3 










6032.64 


13.0 


13.8 
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6.5 


8.1 
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6048.53 
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18.3 
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6054.36 


5.3 








7.7 


5.3 











a The mean wavelength is given. 

b Corrected for the contribution of Aim ,14479.93 estimated from spectral synthesis. 
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